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Abstract 

In this article we analyse, for a simple model, the properties of a practical implementation of a 
fully self-consistent theory of adiabatic large-amplitude collective motion using the local harmonic 
approach. We show how we can deal with contaminations arising from spurious modes, caused by 
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analyse in detail the nature of the collective coordinate in regions where they cross spurious modes 
and mixing is largest. 
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I. INTRODUCTION 



One of the long-standing concerns in nuclear structure theory is to understand how the 
collective properties of a microscopic many-body system emerge from the behaviour of its 
many quantal constituents . Bohr;s liquid drop model can be used to describe many 
collective modes |j| and still provides us with the concepts used in most discussions of 
collective motion. Clearly, a more microscopic method to describe such modes has been 
sought and many methods exist that try to describe collective motion. 

One method commonly used is based on a mean-field approach, the Constrained Hartree- 
Fock-Bogoliubov (CHFB) method [2]. Here a collective energy surface is generated by 
mapping the energy expectation value as a function of the expectation values of a small 
number of one-body operators (mean-field constraints). The principal weakness of CHFB 
is that the collective motion of the nucleus is determined without any dynamical criterion 
for choosing one particular constraint (collective operator) over another. Often the choice is 
inspired by the Bohr Hamiltonian, and one uses only quadrupole and pairing operators. In 
nuclear physics the first applications of the CHFB were to the dominant surface vibrations 
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I which called attention to the need for a more meaningful theory. 
Ideally the method chosen to tackle this problem should allow the dynamics of the system 
to evolve through the microscopic Hamiltonian without the intervention of ad hoc elements 
forcing the system into a specific mode of oscillation. One possibility, followed in this article 
and discussed in detail in the review article '0\, is a method that determines collective 
motion by following the normal modes of oscillation. This incorporates methods such as 
the Quasiparticle Random Phase Approximation (QRPA) Q|. Thus the low-lying physical 
excitations of the system using the small-amplitude harmonic limit (i.e., QRPA) are used 
to self-consistently determine the constraining operator. Restricting to the case of a single 
collective excitation, we then describe a collective path through the energy surface generated 
by the Hamiltonian. 

initial m08t apphcations ignored pacing, and they US ed _ of the Hartree- 
Fock (HF) method. In recent years, three different models including pairing, and thus 
requiring the Hartree-Fock-Bogoliubov (HFB) mean-field approximation, have been investi- 



gated. The first 
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111 ] describes a two crossing levels, with an additional pairing interac- 



tion. One problem first encountered in that work was that of dealing correctly with spurious 
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modes arising from spontaneously broken symmetries within the mean-field approximation 
even though for this model one can solve this problem exactly, unlike the problems discussed 



below. The second 



12 



13j is a microscopic 0(4) model with pairing and quadrupole in- 



teractions which focuses on shape coexistence phenomena in a schematic way. The third 



set of applications [lj, [15( uses the Pairing Plus Quadrupole (PPQ) model for semi-realistic 
nuclei using techniques developed earlier [l6| . These calculations typically use model spaces 
consisting of two major shells for both neutrons and protons, leading to large dimensions for 
the QRPA, and extremely time-consuming calculations. A balance is found by introducing 
a basis of operators which gives a truncated expression of the QRPA in such a way that it 
gives a reasonably accurate approximation of the low-lying excitations This method 

has recently been extended to finite rotational frequencies jis| . 

In this paper we do not make such a truncation, but we will investigate the nature of the 
approximations made before doing the projection on a small basis of operators. To this end 
we employ a single j-shell PPQ model. The main difficulty encountered here is associated 
with admixtures of spurious modes (Nambu-Goldstone modes) arising in the formalism as 
vacuum-degenerate solutions in the excitation spectra which do not correspond to physical 
excitations. Though at extrema on the energy surface spurious modes do not mix with 
other modes, this is not necessarily true at other points on the energy surface, especially 
after additional approximations have been made. Since spurious modes do not correspond 
to physical excitations, any mixing of these modes with other excitations — we will refer 
to this phenomenon as 'spurious admixture' — can potentially lead the system away from 
the collective path. In order to ensure that spurious modes behave correctly far from stable 
equilibrium we need to modify our algorithm slightly, and the study of these modifications 
is the main focus of this article. 

This article is organised as follows: First, in section |H] we give a brief overview of the 
Local Harmonic Approximation (LHA) followed in this article. We also give a description 
of our projection technique to remove spurious admixture from the formalism. Secondl y in 
section ILTT] following the lead given by references Q, Q, Q, Q, Q Q, Q, H, Q Q Hfl , 
the single j-shell PPQ model is presented and solved. We present our main results from 
testing our methods of removal of spurious admixture, as well as discussing a complementary 
investigation of the effect of the pairing strength on the collective potential. In the next 
section, Sec. II VI we then study the behaviour of our techniques at finite angular momentum. 
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Finally, a summary and outlook is given in section 



II. LOCAL HARMONIC APPROXIMATION 

In this section we give a brief overview of the LHA as used in this article; for a more 
detailed exposition of the theory see Ref. We also give a description of two methods to 
remove spurious admixtures from the formalism. 

The LHA starts from assuming a non-relativistic many-body Hamiltonian i? nuc i e ar that 
is capable of describing the low-energy properties of nuclear structure. By means of the 
time-dependent mean field approximation, the quantum problem is turned into a classical 
Hamiltonian problem and the description of collective motion can now be formulated fully 
as a classical decoupling problem. We now try to find collective coordinates such that the 
mixing between collective and non-collective degrees of freedom is small. The Hamiltonian 
contains N — the number of quasi-particle states squared — canonical pairs of coordinates 
£ a G £ = {£\ • • • > and momenta 7r Q G n = {ni, . . . , 7Tj\r} in the Hamiltonian 7i = 
7Y(£,7r). In the adiabatic limit, valid when the collective motion is slow, we Taylor expand 
the Hamiltonian up to second order in momenta. We find 

m,7r)=V(0 + ±7r a B a P(07r p . (1) 

The dynamics of the system are thus characterised by the point function V(^) and by the 
reciprocal mass tensor B a @, which also plays the role of a metric tensor. 

We now change representation from the initial set of coordinates and momenta (£, ir) 
to a new set (q,p) through point transformation, i.e., q = /(C), in the hope that we can 
approximately decouple the collective from non-collective motion. Using the standard Ein- 
stein convention — where a comma denotes a partial derivative / )Q = J^- — the point 
transformation has the form 

r = 9 a (q\---,q N ) = g a (q) , «a = f%P» , (2) 

where g is the inverse of /. The Hamiltonian takes the new form 

H(M=H(q,p) = V(q)+±p li B'»{q)p v . (3) 

Here a bar denotes the transformed potential (V = V(£(q))) and the transformed mass 
tensor (B). 

4 



The basic equation of the problem can now be derived by looking at the dynamical 
fluctuations at an arbitrary point in coordinate space, (the Local Harmonic Equation (LHE)) 



where a semi-colon denotes a covariant derivative. The second equation states that the force 
should be in the direction of one of the eigenvectors, 



For exact decoupling this direction is (co-) tangential to the collective path, but we construct 
an algorithm for finding cases for approximate decoupling by not imposing this as a con- 
dition. Instead, as discussed in Ref. [J, we do have a measure for the quality of collective 
path based on the deviations from this criterion. 

In the equations above, the covariant derivative is defined using the metric B a ^ provided 
by the kinetic energy. The calculation of the covariant derivative requires the inversion 
of the mass matrix, which is fraught with numerical difficulties when B has zero eigenval- 
ues. Therefore we shall make the usual assumption that we can ignore the effects of these 
curvature terms, i.e, V a p ~ V a /3- 

A. Method of removing spurious admixture 

It can be shown that for the HFB theory considered here, the LHE is a simple gener- 
alisation of the QRPA to non-equilibrium states. Thus, as in the usual equilibrium HFB 
plus QRPA framework spurious modes arise as artifacts of mean-field symmetry breaking. 
In a numerically exact calculation, without approximations, such modes decouple from the 
problem, and we so not have to consider them in detail. When approximations are made, 
especially the simplification arising when neglecting the covariant terms in the derivatives, 
the situation is different and we get spurious admixtures. Since these admixtures do not 
correspond to physical excitations we must remove their components completely from the 
formalism without altering the meaning of the collective coordinate. In all cases we know 
the spurious operator; it may correspond to a coordinate or momentum, but we do not 
normally know the conjugate variable. 




(4) 
(5) 



(6) 
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Our method is easily generalised to include any number of spurious coordinates and is 
therefore not specific to the model we solve in this article. As we will see, the main difficulty 
with this approach is that we do not know the conjugate variables and we will therefore 
have to devise a way to approximate an operator for the momenta in order to complete our 
projection method (this is discussed below). 

The first step is thus the construction of a set of approximate conjugate variables. This 
is best illustrated for the particle number, although other operators are dealt with similarly. 

Thus we have q tCe = AA Q , and its conjugate momentum p a = (j) a . We are looking for 
operators that satisfy the conditions 

B a/3 AT^ = , 4> a V a/ 3 oc J\f >a . (7) 

Let us now turn to the canonical coordinates and momenta and the question of how to 
define the momentum variable <p conjugate to the particle number A/". Naive inversion of 
the second equation in (J7j) gives 

<f<xM,a(V,af>)- 1 ■ (8) 

We thus approximate by 




where eP a are components of the eigenvectors and 6j the eigenvalues of the potential matrix. 
Unfortunately V a/ 3 is usually singular due to zero eigenmodes of the QRPA matrix, and we 
need to replace (|8|) by a singular value decomposition. To this end we remove those terms 
in the sum corresponding to eigenvalues that are close to zero €j ~ 0. In practise, the size 
of the excluded zone is found through trial and error. 

Finally, we need to ensure that the operators arising from our construction are all con- 
jugate to each other: In a generalised representation we have a set of canonical conjugate 
coordinates and momenta connected with spurious modes 

q% , Vie {1,2,..., N s } , (10a) 

p? , Wje{l,2,...,N s } , (10b) 

where N s denotes the total number of spurious modes. For well defined coordinates and 
momenta we have q\ a p^ = Sj, but our construction above usually does not satisfy this 
(partially due to neglect of covariant terms, and partially due to numerical problems). 
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We start by defining the overlap matrix 

oj- = M . (ii) 

which, had we have done the computation consistently, would be the identity matrix corre- 
sponding to the Poisson brackets between canonical coordinates and momenta. To achieve 
canonicity we now diagonalise the overlap matrix 

(s-^fojs? = $n , (12) 

where S are a set of orthonormal eigenvectors and 77 are the eigenvalues of the system. This 
allows us to define a new 'tilde' basis 

Ca = . (13) 

P- = ~S}P? , (14) 

T j 

satisfying the orthonormal condition q l a p^ = S l j, used in the definition of the projection 
operator 

pi = €- iJi 

The projection matrix is now applied to the potential and mass matrices, 

Vrt = P£VjpPg , = I* . (16) 

The related QRPA eigenvalue problem has now a number of zero eigenvalues, and the 
spurious modes do not mix with with any other ones. 



(15) 



III. MODEL AND RESULTS 

In this section we solve the single j-shell PPQ model for Adiabatic Large Amplitude Col- 
lective Motion (ALACM) using the methods discussed above to remove spurious admixture. 
The PPQ model we use here is a good starting point despite its limitations as it is relatively 
simple compared to models with more realistic forces and allows us to choose a conveniently 
small basis set with which to work. Furthermore, we would like to 'bridge the gap' between 
simpler applications starting from a partial Hamiltonian and semi-realistic descriptions of 
real atomic nuclei (see discussion above). 
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Following reference 22] the PPQ Hamiltonian will be treated in the Hartree-Bogoliubov 
(HB) framework. We start with the traditional single j-shell PPQ Hamiltonian 



H = - 



(17) 



M 



which is a straightforward two-body interaction Hamiltonian. Here colons denote normal 
ordering and k and \ are the coupling strength of the quadrupole and pairing interactions, 
respectively. We can only impose correct particle number on average, which is done by a 
Lagrange multiplier. When looking at non-zero angular momentum we impose an additional 
constraint for the expectation value of J x . 

The quadrupole and pairing operators are to be written in the simplest form possible in 
the original particle basis 

2 



Q M oc ^ 

mm' 

p = EH 



j—m—M 



J J 

m —m' M 



3 -m, 



(18) 
(19) 



m>0 



where c^ m and c m are creation and destruction operators on the single-particle state \m), and 
the object in square brackets denotes a Clebsch-Gordan coefficient. Since we are working in 
a single shell, we can absorb the reduced matrix element of the quadrupole operator into k, 
and use an equality in Eq. (fTH|) We use the definitions 

'P) /' . (20) 



qM 



Q 



M 



q±i = o, po 



where ang ular brackets denotes the expectation value. We will use the Hill- Wheeler coordi- 



nates 



23] 



(3 = ±Jq 2 + qi 



■ 1 



arctan 



-(12 



(21) 

YVc use the standard HFB formalism as can be found in Ref. 0. Denoting the quasi- 
particle creation operators by b k , we find that for a small fluctuation (labelled by the two- 
quasi-particle index kk') around a given HFB state we have to second order 



H. 



E 

[20] 
l kk' 

•Akk'iv 



H 

bk'bkH 

b k rb k Hb\b\, 

b k 'b k Hbibi> 



(22) 
(23) 
(24) 
(25) 
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Here E denotes the HB vacuum energy, if t 20 ! denotes the gradient of the energy surface at 
that point and the matrices A and B are the usual QRPA matrices, which are simple related 
to the__mass and potential matrices of the LHE. The QRPA equations may be cast in the 
form 



i{A + B)P' i = Q» . 
-i(A-B)Q" = n 2 P» 



(26) 
(27) 



where Q and P are canonical coordinate and momentum variables, respectively. Furthermore 
from equations (J27)) we deduce 



(A + B)(A-B)Q" = n 2 Q^ 



(28) 



which is equivalent to the LHE equation 

As sketched above, there are exact but spurious solutions to the QRPA which are par- 
ticularly related zero eigenmodes of eigenvalue equation (|2~%j). Operators associated with 
rotational invariance in normal space J and in gauge space N, generate new states which 
resemble collective excitations, and thus 



H, J x 



, 



H,N 

H, Jy 



, 



-uiJ 7 



H,J Z 



uiJo, 



(29) 



The last two terms are zero for the J = ground state, but at finite u>, iJ y and J z form a 
canonical pair, a coordinate and a momentum. In that case the set of spurious coordinates 

i S r' = {Af,j x ,j,}. 



A. Results 



We have first investigated the quality of the projection scheme for the single j-shell PPQ 
model. First consider axially symmetric (q±2 = 0) ALACM, which corresponds to following 
the lowest QRPA mode as a function of the collective coordinate Q. The model space is 
chosen to be small: j = | and N = 4 fermions and with interaction strengths K — 1 and 
X = 0.333. We have calculated the mean-field parameters along the collective path. The 
results are shown in figure H 
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FIG. 1: ALACM following the lowest QRPA mode as a function of the collective coordinate Q. 
Figure^ shows the collective potential V, figure ^ shows the /3-deformation and figured shows 
the pairing parameter p$. Each panel shows the naive algorithm (dotted line) and the zero-mode 
corrected one (solid line). The scale of all displayed quantities is arbitrary, for discussion of these 
see main text. 

For the corrected algorithm (solid line) we identify two minima and three maxima at 
Q ~ {—2,0} and Q ~ {—5,-1.25,2}, respectively. The lowest-energy solution at the 
starting point Q = corresponds to a prolate minimum and the second minimum at Q ~ —2 
is oblate. Comparison of the collective potential with the deformation parameter (3 and 
pairing parameter p exhibits the well-known competition between pairing and quadrupole 
forces. At the boundaries of the collective coordinate Q ~ { — 5,2} the system shows maximal 
oblate and prolate deformation, respectively, but the pairing field collapses, po ~ 0. The 
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FIG. 2: The overlap of the collective operator with the quadrupole, pairing, and particle-number 
operators as a function of the collective coordinate Q. The solid line shows Xq, the dotted line 
Xp, and the dashed line X/v, Eq. (|3U|) , Figure^ displays the results for the uncorrected algorithm 
and figure \Bp for the corrected algorithm. 

converse is true at the central maximum where the system has collapsed into a state of 
pure-pairing; q±M ~ 0. 

Figure H shows clearly that without correctly removing spurious modes (dotted line) the 
ALACM algorithm fails outside a limited range of the collective coordinate —3.75 < Q < 0.5; 
we note that the collective potential and associated mean-field parameters follow a different 
path from the corrected algorithm (solid line). To better understand this we turn to an 
analysis of the collective coordinate and calculate the degree of overlap of the collective 
operator with the operators contained in the Hamiltonian: O G {Qo, (P^ + P)i X}, 

X = J2^0 a . (30) 

a 

The results are shown in figure El 

For the naive algorithm shown in figure |2k (c./. figure (TJ;) we see that the collective 
coordinate is dominated in the region of maximal pairing by a quadrupole- like operator 
[T ~ Qo). In particular, where Q ~ —1.25 the pairing-like parameter Xp tends to zero, 
and we have a pure j3- vibration. From figure |2t we see clearly why the naive algorithm 
fails at Q ~ —1.5; first, we see a sudden increase in the parameter Xn, and secondly, a 
sharp change in the pairing-like parameter Xp. From this it is evident that the the spurious 
coordinate associated with particle number M is no longer sufficiently decoupled from the 
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FIG. 3: QRPPA energy eigenvalues squared (dotted lines). The solid line shows the QRPA mode 
selected by the path following algorithm. Figures OK and Ot display results for the uncorrected 
algorithm and figures Ob and OJi display corresponding results for the corrected algorithm. The 
scale of all displayed quantities is arbitrary, and determined by the value of k = 1. 

collective coordinate. 

This can be seen in figure 01 where we study Q 2 , Eq. (|28|1. On the left we show the results 
of the naive algorithm, and on the right the results from the zero-mode corrected one; the 
upper two figures differ from the lower two only in scale. It can be seen from figure OK that at 
Q ~ —1.5 the energy of the selected eigenmode changes abruptly. At this point we note that 
the eigenvalue of the collective mode Q 2 and that of the mode associated with the particle 
number operator Q 2 N have become degenerate and an avoided crossing of the two levels has 
occurred. Here the path following algorithm selects the wrong eigenmode, i.e. a spurious 
mode rather than the collective mode, leading the system away from the collective path - 
this is seen most clearly in figure EJ In figure Ob is shown the corrected result whereby the 
collective coordinate smoothly crosses the level associated with particle number. 

In order to illustrate the source of the problems in the naive algorithm we study both 
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FIG. 4: Analysis of the eigenvalues of the QRPA. Top boxes shows eigenvalues of the QRPA (dotted 
lines) following the lowest QRPA mode (solid line) as a function of the collective coordinate Q. 
Diamond symbols indicate computed points (joined by linear interpolation). Bottom box shows 
the measure of the overlaps Xo- 

the real and avoided crossings present in figure Ek- We see that the avoided crossing has a 
large admixture of particle number, and that even the real crossing has a smallish admixture 
(indicating it is probably a very narrow avoided crossing). We have further investigated the 
nature of the avoided crossing by using a variety of step lengths. We find that the result is 
not influenced by changing the step size, and that we can not take steps so large that we 
run straight through the avoided crossing. 

We have investigated changes in the collective potential with respect to variations of the 
pairing strength \ from 0.28 to 0.34. The results are shown in figured From figured we see 
that the collective potential is higher in energy for lower pairing strength and with a higher 
maximum relative to the ground-state. Both ends of the collective potential terminate at the 
same energy, since each of the five systems are identical, that is, where only the quadrupole 
moment is in effect yielding maximally deformed axially symmetric spheroids. Moreover, 
we observe from figures Eb-Et that the maximum pairing and /3-deformation attained by the 
system is identical in all cases of pairing strength x as is expected. 
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FIG. 5: The result of varying the pairing strength. See figure^for an explanation of the subgraphs. 

We conclude that the zero-mode-corrected algorithm seems to be able to deal much more 
efficiently with crossing modes than the naive one. The results obtained with this algorithm 
have all the properties of the full results, and do not suffer from spurious admixtures and 
avoided crossings. Since finally we want to understand what happens for rotating nuclei, we 
must look further at collective motion at non-zero angular momentum. 

IV. PPQ MODEL AT FINITE ROTATIONAL VELOCITY 

To perform a representative test-case of ALACM at finite rotational velocity we have 
minimised the mean-field energy, and computed the QRPA eigenvalues for values of angular 
momentum < J x < 6, where J x = ^J^. The results are shown in figure |H1 

Figure |H1 shows the behaviour of the ground state and its QRPA excitations = Eo + fl^, 
where E is the ground-state energy and is the /ith QRPA frequency. As the angular 
momentum increases, the pairing field weakens and collapses to zero for J x = 3. Until that 
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-2.5, 



1 



c 

— V^co 

— E 



2 3 4 5 
<J > 



FIG. 6: Partial energy spectra of the PPQ model at finite rotational velocity. Figure shows 
ground-state energy and a selection of QRPA eigenvalues for fixed quadrupole strength k = 10, 
fixed pairing strength, x = 0.333, and rotational velocity < J x < 6. The scale of all displayed 
quantities is arbitrary. 



collapse the ground-state (solid line) is degenerate with two eigenmodes associated with 
the two spurious operators M and J x . Thereafter, the only spurious operator is J x . The 
dashed line corresponds to a solution of the QRPA JF W that appears at the value of the 
rotational frequency squared (u^), corresponding to the conjugate pair of operators J y and 
J z . Most interesting here, is the appearance of the lowest-lying excitation corresponding 
to a collective mode T c (dot-dashed line). As the rotational energy increases the collective 
excitation comes closer to the ground-state. 

Until pairing collapse modes mix between pairing and deformation modes. At the point 
of collapse such modes decouple and we should really solve the ordinary RPA. Nevertheless, 
even solving QRPA we can easily distinguish between the excitations that change particle 
number and the deformational modes (continuing lines) which continue to be well behaved 
into the realm of no pairing. These modes correspond to excitation levels of the standard 
Nilsson model 
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A. Model results at finite rotational velocity 

We study angular momentum J x = 2, such that the ground-state of the system has non- 
zero pairing and quadrupole effects (c.f. figure EJ), we have computed ALACM for the single 
j-shell PPQ model at finite rotational velocity following the lowest collective mode of the 
QRPA. For illustrative purposes we will in the first instance remove only admixtures with 
the particle number operator. 

The quantities Xo G {Xq, Xp, Xn, Xj x , Xj z } give the measures of overlap for the col- 
lective operator with the quadrupole, pairing, particle number and rotational operators, 
respectively. 

The results are shown in figure [7| As can be seen, the computation has produced values 
for —0.6 < Q < 0.3 only, due to an instability arising in the path following algorithm 
discussed further below. For negative values of the collective coordinate we see a decrease 
in /3-deformation (J2b) and an increase in 7-deformation JZt), indicating that the system 
is becoming more triaxial, while the pairing parameter po increases as well (01). We thus 
expect the components of the M = ±2 quadrupole operator to be important. 

As can be seen from figure [7^ a spurious admixture connected with J x is causing problems 
for Q < 0. As can be seen from figure [7f the eigenvalue of the collective coordinate Q c (solid 
line) and the eigenvalue of the coordinate associated with the spurious mode (dashed 
line) are approximately degenerate. There are therefore two approximately degenerate paths 
for the path-following algorithm to distinguish, one of which is a spurious solution. 

We note one interesting factor outstanding from our discussion of the results thus far. It 
can be seen that the collective coordinate represented in figure[7f runs approximately parallel 
with the next higher lying state and though we are unsure of the physical significance of 
this, it would be interesting to analysis the collective and non-collective coordinates in more 
detail. To this end, better technology aimed at giving a more detailed description of the 
system will be introduced below. 

We now project out all components associated with spurious modes M , J X) and J Z) and 
secondly, we will compute the measures of overlap associated with all operators of the PPQ 
model. 

We work with the four quantities O G {Q , (Q+1 + Q-1), (Q+2 + Q-2), (V*+V)}. Hereon, 
X e {X Jx , X Jz ,X N } — the quantities connected with projection — are zero eigenvalues of 
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FIG. 7: ALACM for the PPQ model at finite rotational velocity following the lowest QRPA mode 
as a function of the collective coordinate Q. The scale of all displayed quantities is arbitrary. 

the QRPA eigenvalue equation and otherwise will not be alluded to. The results are shown 
in figures IHHni 

From figure |H1 we see that for positive values of the collective coordinate the collective 
path has been successfully computed to its end. This is identified by pairing collapse po ~ 
for Q ~ 0.8 (figure IBtO at a maximum in the collective potential (figure [H^). As the system 
approaches pairing collapse the /3-deformation approaches its maximum value and the 7- 
deformation approaches a minimal value (figure IHb). In other words, at pairing collapse the 
system is maximally axially deformed with some small admixture of triaxial deformation. 
For < Q < 0.8 the collective coordinate becomes increasingly more pairing-like for greater 
values of Q and less quadrupole-like. In particular, this statement is true of all components 
of the quadrupole operator, i.e. for Xq , Xq 1 and Xq 2 , (solid, dotted and dashed lines 
respectively) and thus the collective coordinate becomes very pairing-like trying to force the 
system back toward equilibrium. 

For negative values of the collective coordinate the results show an instability at Q ~ 
—1.25, where the computation has failed. From figure|Hh we see that the collective coordinate 
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FIG. 8: ALACM for the PPQ model at finite rotational velocity following the lowest QRPA mode 
as a function of the collective coordinate. The scale of all displayed quantities is arbitrary. 

Q c (solid line), becomes degenerate with a second coordinate, labelled f2 6 (dashed line). 
Thus the assumption of a single collective coordinate fails here. This is quite common in 
this approach, and requires and extended algorithm. 

In order to convince ourselves that this is the correct physics, and we are not looking at 
the spurious admixture problem again, we analyse the QRPA modes f2 M in some detail. We 
have computed the measures of overlap for a few of these modes. At equilibrium there are 
three spurious modes, and the fourth eigenmode (next highest in energy) is the collective 
eigenmode we are following Q c . We concentrate our efforts on the next few QRPA modes, 
namely, the fifth Q 5 (dashed line), sixth Q e (dot-dashed line), and the seventh Q 7 (dot-dot- 
dashed line, in figure |Hli). We note from figure |HF that these modes do not remain in this 
order for all values of the collective coordinate since there are crossings of energy levels in 
the energy spectra. The labelling of these modes however, shall remain fixed as defined at 
Q = 0, regardless of their subsequent order. The results are shown in figure |H1 

We have computed the measure of overlap 

*o = X^°>* ' V/iG {5,6,7} . (31) 
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FIG. 9: Analysis of low- lying eigenmodes of the QRPA eigenvalues of the QRPA following the 
lowest QRPA mode as a function of the collective coordinate. The scale of all displayed quantities 
is arbitrary. (For discussion of displayed quantities see main text.) 

Figure El shows the results of our analysis for Q 5 (figure Ek), ^6 (figure Eb), and Q7 (figure 
Et). In each sub-figure is shown the measures of overlap Xq (solid line), X Ql (dotted line), 
Xq 2 (dashed line), and X P (dot-dashed line). 

The fifth eigenvalue (figure Eb) and seventh eigenvalue (figure Et) both exhibit a complex 
arrangement of all measures of overlap, whose evolutions are not especially discernible; we 
will return to this later. The sixth eigenmode, figure Eb, contains exclusively admixtures of 
Xq 1 and Xq 2 meaning that this is a mode that drives both 7-deformation (Q±2) and tilting 
(Q±i). At the point of instability, Q ~ —1.25, the tilting character is particularly strong. 



V. SUMMARY AND OUTLOOK 



In this article we have investigated ALACM using the LHA for the single j-shell PPQ 
model. We have computed the collective path, as well as a variety of characteristics of the 
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mean-field as a function of a single collective coordinate. In order to analyse the collective 
coordinate in more detail we introduced the measures of overlap, which turned out to be 
an indispensable tool to analyse the collective coordinate especially in regions where en- 
ergy levels become degenerate. We introduced a projection technique to remove spurious 
admixture, which works by projecting a subspace out of the computed mass matrix and 
potential matrix free from spurious admixture. This method is easily extended to remove 
the admixture of any finite number of spurious modes. The mixing of spurious modes with 
the collective mode brings to light an inadequacy of our initial assumption that all modes 
are approximately decoupled from the collective coordinate. 

Using the LHA, collective motion was shown to be non-linear in nature and dominated by 
neither pairing nor quadrupole degrees of freedom exclusively. Rather, the system self-selects 
either pairing, or quadrupole, or an admixture of the two to dominate the dynamics and 
thus select how best the system deforms. More especially, we have found the remarkable 
feature of a shorter collective path for stronger pairing strength which is not an obvious 
result. 

This leads us to an important limitation, namely, that a single collective coordinate is 
not sufficient to describe ALACM in regions where the collective coordinate undergoes real 
crossings with other physical solutions which mix with the collective coordinate. Further- 
more, in the course of this investigation we have stumbled across the apparent importance 
of tilting solutions arising in the quadrupole moment. 

We have not computed the covariant derivative of the potential, but we have rather com- 
puted a truncation of this which leaves a normal derivative and losing some terms. These 
terms may well be more important than we suspected at the start. This is something we 
would like to investigate further. In addition to this, we have attempted to work within 
the full configuration space available in order to allow a large set of degrees of freedom 
to contribute to the dynamics. It is worth mentioning here that work on a semi-realistic 
approach to ALACM using the PPQ model and incorporating many j-shells uses a trun- 
cated QRPA which inadvertently avoids handling spurious modes directly. The projected 
QRPA was introduced as a computational aid where it was assumed that coupling of normal 
and spurious modes with the collective coordinate would be small. As we have shown, this 
assumption may not be correct. 

Finally, we note that the LHA is not restricted to a solution for a collective path only, 
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i. e. we could take into account more than one collective coordinate to describe the dynamics 
of ALACM. The next step is to modify the algorithms we have built, such that they are 
more apt to describe nuclei at finite rotational velocity using two collective coordinates to 
determine a collective surface. 
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